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A NEW METHOD OF SOLVING NUMERICAL EQUATIONS. 



BY DK. H. EGK3ERS, MILWAUKEE, WISCONSIN. 

Among the methods for solving numerical equations, there are such as 
teach, to find from one or more initial approximate values a more approxi- 
mate one of the root required. Such methods are, e. g., the method of New- 
ton and the rule of false position. These methods become in certain cases 
impracticable because of slow convergence. The Newtonion method, e. g., 
converges very slowly toward a multiple root. For in the case of a sim- 
ple root the convergence is of quadratic order, while in the case of a multi- 
ple root the acceleration is only linear. That is to say : in case of a simple 
root every correction is proportional to the square of the preceding cor- 
rection ; but in case of a multiple root every correction is proportional to 
the first power of the preceding correction. 

It is true, the multiple roots can always be separated from the given equa- 
tion, but the troublesome case of nearly equal roots remains. Besides this, 
the labor of separating multiple roots is not always small, and it is there- 
fore desirable to avoid it, if possible. 

The object of this essay is, to communicate a method of calculation 
which is independent of the relations of roots to one another, to wit, a 
method which converges with the same rapidity whether the roots be sim- 
ple or multiple: The acceleration towards the roots is in every case of 
quadratic order. This algorithm is the following : 

Theorem. Liekflx) = be a numerical equation of the n th degree, x an 
arbitrary approximate value of one of its roots, f(x) and f'{x) the first and 
second differential quotients, then the next corrected value x x is found to be 

(1) x —x /fooXffco) 
or briefly 

( 2 ) x 1 =x — (jyl_ff> = -^o)- 

The next corrected value x 2 would be the same function of x x as x x is of 
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Now let x be one of the roots of the given equation f(x) = 0. In this 
case we would derive from equation (2) x x =x; for the correction would 
vanish, as f(x) is a factor of it. Further, let us suppose x to be an ap- 
proximate value of this root, so that 

Xq = x -J- e, 
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then the corrected value x t would no more be equal to x , but we would 
have, according to (2): 

«! = *(* + «) = F(x) + e.F(x) + ieW'ix) + l<?F"\x) + 

and x x would differ from the root by some quantity e x , therefore 

x 1 =x + e 1 = F(x) + e.F{x) + le 2 F"{x) + 

Now x = F(x), consequently we get an expression for e x , 

e x = e.F(x) + %(?F'(x) + \#F"{x) + 

The function F(x), the first differential quotient of F(x) possesses the 
property (and this is essential for the establishment of our method) that, the 
first differential quotient, F(x), vanishes for f(x) = 0, even in the case when 
the root x is a multiple root of any order. Therefore the error e x of the cor- 
rected value of the root x x is 

e x = ie*F'(x)+^F"(x)+ . . .; 
to wit, the new error e x is proportional to the square of the preceding error, 
whether the root in question is simple or multiple; that is, the acceleration 
of our method is always of second order, whereas the method of Newton, 

in case of a multiple root furnishes a convergence of linear order. 

If it should happen that F"(x) vanishes for f{x) = 0, the acceleration 
would be of the third order. 

This algorithm is only one particular case of an infinite number of meth- 
ods of every possible rate of convergence. But the amount of work con- 
nected with those algorithms of greater convergence is not fully balanced 
by the greater rapidity of convergence. For example, an algorithm of cubic 
acceleration, for simple or multiple roots, is the following : 

x -x - f(f 2 -ff") 

where/,/,/' are functions of x ,f the first, and /' the second differential 
quotients of the given equation f(x) = 0. 

As it would require considerable space to prove the above statements, I 
suppress the demonstration here, but I may be permitted to add one exam- 
ple: 

Ex. Given the equation (x — 2) 3 (x — 4) = with the multiple root 2. 

Taking for initial value x = 3, we find successively 

x x = 2.5, x 2 = 2.07, x 3 = 2.0009, x A = 2.0000001, etc. 

Observation. Always that root is found the modulus of which is the near- 
est to the modulus of the arbitrary initial trial value. If in the foregoing 
example we would commence with x = 5, we would approach the root 4. 



